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INTRODUCTION 

Much theoretical and experimental attention has centered on the photo- 
chemistry and energy balance of the plasma in the lower thermosphere [Hanson, 
1963; Norton et al., 1963; Dalgamo ot al., 1963; Holmes et al., 1965; Spencer 
et al., 1965; Swider, 1965; Donahue, 1966; Bauer, 1966; Johnson et al., 1967; 

Brace et al., 1969; Brinton et al., 1969] . A shortcoming of most studies is that 
the investigators considered the problems separately, using ion-composition 
profiles obtained at one time along with electron concentration and temperature 
profiles taken at a significantly different time or location, assuming models for 
iinmeasured parameters. Because many photochemical and thermal processes 
in the atmosphere are interdependent, near-simultaneous measurements are de- 
sirable in order to produce a consistency of results when studying the physics 
of the thermosphere. 

This paper presents a set of F region vertical profiles of positive ion com- 
position and electron concentration and temperature obtained above Wallops Island, 
Virginia, on August 26, 1966. These meas\irements resulted from experiments 
flown on Nike- Tomahawk rocket NASA 13.06 at 1351 EST (1851 UT). Measure- 
ments of molecular nitrogen concentration and temperature [Spencer et al., 1969] 
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(Thermosphere Probe NASA 18.05, 1331 EST) and spectroscopic measurements 
of a few lines in the extreme ultraviolet range [ Hinteregger, personal communi- 
cation] (NASA 4.101, 1313 EST) on the same day just before iuC NA.SA 18.06 
flight complete the set of data used to study the photochemistry and plasma 
energy balance simultaneously. This analysis considers the production and 
loss of ions and electrons in the F-region, as well as the spectrum of ionizing 
radiation, ionization cross sections, plasma cooling processes, neutral and 
charged r<article concentrations, temperatures, and recombination coefficients. 
The results include rates of the primary chemical reactions responsible for the 
formation of the ionosphere as well as those for the electron and ion cooling. 

INSTRUMENTATION 

The Ion Spectrometer 

As the operation of Bennett mass spectrometers has been described by 
Taylor et al. [1963, 1965], this paper contains only details of this particular 
experiment. The sensor, a 7-5 cycle cylindrical ceramic tube with 3 mm grid 
spacing, operated at a frequency of 5.4 MHz. The spectrometer tube was 
moimted coaxially along with the associated electronics and earth telescope (a 
part of the optical -aspect system), in one end of the thermosphere probe experi- 
ment package. The ion mass range from 12 to 36 amu was scanned every 0.15 
seconds. 
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The Cylindrical Electrostatic Probes 


» 


Spencer et al. [1965] have described the derivation of electron temperature 
(T^ ) and concentration (N^) from electrostatic probe measurements. Two 

a 

probes made independent measurements on this flight: the probe sensors had 
identical dimensions (collector: diameter = 0.022 inches, length = 9.0 inches; 
guard: diameter = 0.067 inches, length = 5.0 inches), but each was made of a 
different material to permit a study of the electrode surface/plasma interactions 
upon the measurements. One was made of gold-plated stainless steel; the other • 
was also gold-plated, but made of beryllium copper. . 

RESULTS 

Electrostatic Probe Data 

The volt-ampere characteristics from the two probes agreed well through- * 
out the flight. However, the data showed evidence of local contamination, per- ^ 
haps caused by rocket motor outgassing, during approximately the first 30 
seconds after the instrument package was ejected from the nose cone. Con- 
sequently, the results below 150 km came from the descent data which did not 
exhibit this problem. Figures 1 and 2 show the altitude profiles of and T^ 
respectively. Figure 2 also shows the molecidar nitrogen temperature (T^ ) 
profile, derived from the Nj scale heights measured on NASA 18.05, just 20 min- 
utes earlier, illustrating the difference in T^ and T^. 
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Ion Composition 


The total Ion current detected during each spectral sweep was normalized 
to a vertical profile of total ion concentration (N^ ) to derive the concentrations 
of the individual ion constituents. This profile was established assuming 
N. = N, and by normalizing (at F-maximum) the smoothed profile of electron 
concentration (N^) measured by the electrostatic probe to a curve resulting from 
groimd based-lonosonde measurements from 130 to 240 kilometers. The electro- 
static probe thus provided an extrapolation of the ionosonde profile from F-maximum 
to peak altitude. Before normalization to the total ion concentration profile, the 
individual ion-current measurements were corrected for mass discrimination, 
using factors derived from laboratory data. After the normalization process, 
we eliminated the wake data; the respective ascent and descent data then agreed 
within 20 percent. Averaging these remaining points produced the final profiles 
shown in Figure 1. 

DISCUSSION 

The Electron and Ion Thermal Balance 

By employing the profile of T^ measured by the electrostatic probe experi- 
ment on this fli^t to calculate the electron cooling rates as a function of altitude, 
the derived electron and ion thermal balance can be compared to earlier results 
of measurements by Spencer et al. [1965] and Brace et al. [1969] . In the 
derivation the procedure described by Brace et al. [1969] was used except that 
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the effects of electron heat conduction (following Banks [1966]) were also in- 
cluded. The ion temperature (T^ ) was calculated by equating the electron-ion 
heating rate to the ion-neutral cooling rate. As mentioned earlier, Nj tem- 
peratures (see Figure 2) and concentrations were obtained from measurements 
reported by l^)encer et al. [1967, 1969] . The Jacchia 1964 model atmosphere 
provided the remaining neutral constituent concentrations, using the measured 
neutral temperature to determine the appropriate model. Figure 3 shows the 
calculated local electron cooling rate (solid line) and the local heating rate 
(dashed line). The difference indicates the amount of heat conducted into or out 
of the region, showing that above 210 kilometers heat conduction acts as a heat 
source, below a heat sink. The net effect is that the F^ region peak of local 
cooling shifts to a higher altitude than the peak of local heating. At altitudes 
below 160 kilometers the rates increase toward higher values, and exhibit 
maximum values in the E-region. The cooling rates, throughout the entire range 
of altitude, resemble those derived from data obtained at Wallops Island in 1960 
and 1961 [Spencer et al., 1965] at similar levels of solar activity (F^^ , - 125 
vs Fjg 7 = 103 for the 18.06 flight). 

Ion Chemistry 

The basis of our photochemical analysis was the ultraviolet spectrum and 
the absorption and ionization cross sections given by Hinteregger et al. [1965] . 

O 

The few lines in the EUV spectrum monitored above 300 A only 40 minutes 
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before this flight [Hlnteregger, personal communication] did not exhibit any 
consistant trend in intensity variations. However, at the time of the charged 
particle experiments, the 10.7 cm flux (which reflects die level of solar activity) 
exceeded by 50 percent its value in 1963 when Hlnteregger's measurements 
were made. Considering that the low wavelengths in the EUV spectrum are 
especially sensitive to variations in solar activity, we scaled the photon fluxes 

0 o 

linearly from a factor of one at 300 A to a factor of three at 100 A. This pro- 
cedure, which keeps the spectral intensities within experimental uncertainties 
of the earlier EUV measurements, resulted in good agreement with ion com- 
position measurements. We assumed that each 25 ev of photoelectron energy 
produced one electron/ion pair, giving an efficiency factor of 


Eff 



25 


where E^^ is the energy of a photon dn ev) for wavelengths smaller than 350 A; 
0 

above 350 Awe assumed that the efficiency factor was unity. Figure 4 shows 
ionization rates for O^, and (including primary and secondary ioniza- 
tion) calculated by using Jacchia's 1964 model for the unmeasured neutral 
composition and the measured exospheric temperature (950”K). 
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N/ Iona 


Figure 4 illustrates the production of Ions photoionization of N,. The 
loss of ions results primarily from the reactions 


Nj* + 0 - NO* + N 


N/ + Oj - Oj* + Nj 
[Ferguson, 1967], and 

Nj* + e - N + N 



a) 

( 2 ) 

(3) 


The last reaction is electron-temperature (T, ) dependent, (Biondi, 1967] with 


k 

N 


♦ 

2 


e 



where 


k . = 2.5 X 10“^ cm* sec"* 

«2 -0 

Figure 2 shows the T^ altitude distribution. The N,* continuity equation, using 
the photoproduction rates of N,*, the measured concentration of N^, NO*, and 
O * and reactions 1-3, was solved Inr varying k and k . to provide an 

N, 0 N,0, 
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optimum fit to t!i<3 obsorved profiles. The resulting values are 


k . = 8.0 X 10“** cm* sec"* 

Nj 0 

and 


k 

N 


2^2 


6.2 X 10“ * * cm* sec 


These values are lower than those summarized Ferguson [1S07] by about a 
factor of two, which is approximately the uncertainty for the ionization cross 
sections [ Hinteregger, personal communication] . 


NO* Ions 

Production of NO^ ions is by reaction (1) and by reaction 

0* + N - NO* + N r*« ♦ 1 

L o "2J 

and 

N* + O - NO* + 0 fk ♦ 1 (5) 

L " 


[Ferguson, 1967]. Loss of NO* is by dissociative recombination 


with an electron-temperature coefficient of 


where 



C7) 


k = 7 X 10“’ cm^sec"’ 

NO%o 


[Biondi, 1967]. 

Because the concentration [ Jacchia, 1964] is at least a factor of five lower 
than the measured concentration and the measured c<mcentration is nearly 
two orders of magnitude lower than the measured concentration over the en- 
tire altitude range (Figure 1), we can neglect reaction (5). With the knowledge 

of the previously determined coefficient, k ^ , we used (1) and (4) and our 

Nj o 

measurements to evaluate the following coefficient for reaction (4) 

k . = 6.5 X 10"^^ cm^sec"’ . 

O Nj 

Laboratory measurements of this coefficient gave an uf^r limit of 6.7 x 10“ ’ ’ 
cm^ sec~’ [Talrose et al., 1962] . Measurements by Fehsenfeld et al [1965b] 
resulted in a value of 3.0 x 10" cm^ sec“ Schmeldcopf et al. [1966] re- 
ported that this rate may sharply depend on the vibrational temperature of N 2 , 
and this may explain some of the higher values reported for k ^ . 

O Nj 
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O* Ions 


The ions are produced photoionization and by the charge exchange 
reaction 

. ( 8 ) 

Ferguson et al. [1965] estimate the reaction coefficient for (8) to be smaller 
than 10"^ ^ cm^ sec" Considering the density when compared to the O* 
density (Figure 1) that is involved in the subsequent loss reaction, we may neglect 
( 8 ). 

The ions are removed by reaction (4) and by 

. 1 ( 9 ) 

°2j 

[Ferguson, 1967]. 

Using the production rate for O* (Figure 4) and the previously derived rate 
coefficient, k + , we found that a coefficient for reaction (9) of 

O Nj 

k . = 5.9 X 10"^^ cm^ sec" ^ 

0 Oj 

agrees best with the measurements. For this rate coefficient the values 
that occur in the literature vary between 1.8 x 10"^^ [Langstroth and Hasted, 

1962] and 4 x 10" ^ ^ cm^ sec" ‘ [Fehsenfeld et al., 1965a] . 
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Oj Ions 


O, ions are produced pbotolonization (rate sho^vn In Figure 4), by reac- 
tions (2) and (9), and (according to Fehsenfeld et al. [1965c]) by 


+ Oj - 0 * + N 

Removal of is dissociative recombination 

o; + e - 0, 



with 


k 

o 


+ 

2 


e 



where 


ao) 


ai) 




= 3.0 X 10“’ cm^ sec“^ 


[Biondiy 1967] . Using the reaction rates k ^ and k ^ derived earlier, we 


0 o, 


Nj Oj 


can estimate the coefficient for k ^ , assuming photochemical equilibrium. 


N O, 


Our analysis has shown that its value should not be larger than 


k . = 3 X 10“^® cm* sec"* 

N Oj 
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This is lower than the laboratory values summarized by Ferguson [1967], which 
vary from 4.6 x 10" ^ ° to 6.0 x 10" ^ ® om^ sec" ^ 

The production mechanisms (2), (9), and (10) and the ionization rate for 
decrease with altitude, whereas the reoombinatian rate (which is proportional 
to ) tends to increase. From this, although we expected a decrease of 
above 130 km, we found an essentially constant concentration below 160 km. 

A possible e3(planation is that the photoionization rate for is actually lower 
than we calculated. Then the total production rate would decrease significantly 
at lower altitudes where the other production mechanisms involving O^, N^, and 
(which decrease with altitude) become less significant. 

In order to check the consistency of the derived reaction-rate coefficients 
we solved the coupled continuity equations for O^, NO^, and and the 
previously discussed ion chemistry- reactions for the ion altitude distributions 
assuming photochemical equilibrium and using the derived reaction rate coef- , 
ficients. The measured ion profiles diverged from the calculated distributions 
above 230 km, indicating that diffusion processes are important above this 
altitude. We did not calculate the because this procedure Lavolves the He^ con 
centration, which was not measured in this experiment. However, the measured 
N* was used where applicable. 

Figure 5 shows good agreement between measured and calculated ion dis- 
tributions below 230 km. The only significant exception is at lower altitudes, 
which probably indicates that the photoionization cross sections adopted for O* 
were too high. 
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Distributions of NO * and Oi,^ 


Figure 5 shows the measured distributions of NO^ and to be nearly 
parallel; in near agreement with the calculated model. These measurements 
confirm, to an altitude of 240 km, the observation of Holmes et al. [1965] that 
when NO^ and are in chemical equilibrium, the scale heights correspond to 
the scale heights of N, and O, respectively. Above 240 km a sharp increase in 
the NO ^ and O,*^ scale heights is evident. This increase is considered too great 
to be consistent with the rise in gas temperature in this region. Brinton et al. 

[1969] , suggested that diffusion (and consequently plasma temperatures and 
polarization fields) governs ion distributions at higher altitudes. It is important 
to consider the effects of the variability of the neutral composition, the plasma 
temperatures, and the UV spectrum of the sim upon the ion composition. To 
investigate this we compared the characteristic times for diffusion and the 
chemical reactions involving these ion species. The diffusion time is 


T — 1Q-7N (12) 

D 

where H is the scale height of NO'*’ or , and the diffusion coefficient, D, is 
given by 


■ a3) 

en 
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where 9 Is the drag coefficient cgs), and N is the total atmospheric 

dmisity. The production time for NO^ is 


T = 


[NO*] 

k„. (o‘) IH^ 


7 X 10' (sec) , 


a4) 


and the recombination time is 


T=- ~ 3x10' (sec). (15) 

*S<0+eNe 

The altitude where diffusion and chemical reactions are comparable is about 
270 km where N = 7 x 10^ (cm~^). Above this height diffusion should dominate the 
NO* distribution. 

In diffusive equilibrium, assuming for simplicity that T^ is constant with 
height, the diffusion equation for NO* is 


1 3 [NO*] 

[NO*] 3^ ° ■ k T, ■ T, N. 3i 

where 01 ^^+ is the mass of NO*, and Tj is the ion temperature. From (12) it 
follows that in the region where diffusive equilibrium prevails (above the F, 
peak) the scale height will increase because of the polarization field (second 
term in 16). Hence the observed change in the NO* scale height and the O* 
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scale height using similar arguments seems reasonably well understood In 
terms of the transition from chemical to diffusive processes. 
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FIGURE CAPTIONS 


Figure 1. Altitude profiles of measured electron and ion concentrations. 

The ion concentrations are an average of ascent and descent data, derived by 
normalizing the ion currents to the total ion concentration, , which is derived 
from ground based ionosonde measurements and above 240 kilometers , the 
electrostatic probe data. Below 150 kilometers only the descent electrostatic 
probe data are shown. 

Figure 2. Altitude profiles of and derived from experimental measure- 
ments on NASA 18.06 and NASA 18.05 respectively. 

Figure 3. Altitude profiles of derived local electron cooling and heating rates. 
The difference indicates the heat conduction into or out of the region. 

Figure 4. Calculated photoproduction rates for , and 0*, 

Figure 5. Comparison of altitude profiles of ion concentrations calculated 
assuming photochemical equilibrium and derived reaction rates with the ion 
concentrations measured by the ion spectrometer aa NASA 18.06. 
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Figure 1-Altitude profiles of measured electron and ion concentrations. The ion concentrations 
are on averoge of ascent and descent data, derived by normalizing the ion currents to the total 
ion concentration, N., which is derived from ground based ionosonde measurements and above 
240 kilometers, the electrostatic probe data. Below 150 kilometers only the descent electro- 
static probe data are shown. 
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Figure 2-Altitude profiles of Te and Tg derived from experimental measurements 
on NASA 18-06 and NASA 18-05 respectively. 
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Figure 4-Colculated photoproduction rotes for 




KILOMET 


200 


□•AO CALCULATED 

ION SPECTROMETER 

MEASUREMENTS 


ION CONCENTRATION (lONS/CM^) 

Figure 5-Comparison of altitude profiles of ion concentrations calculated assuming photo* 
chemical equilibrium and derived rejction rates with the ion concentrations measured by the 


spectrometer on NASA 18-06* 






